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Abstract 

We review our recent experimental results on the electronic nematic phase in electron- and hole-doped 
BaFe 2 As 2 and FeSe. The nematic susceptibility is extracted from shear-modulus data (obtained using a 
three-point-bending method in a capacitance dilatometer) using Landau theory and is compared to the 
nematic susceptibility obtained from elastoresistivity and Raman data. FeSe is particularly interesting in 
this context, because of a large nematic, i.e., a structurally distorted but paramagnetic, region in its phase 
diagram. Scaling of the nematic susceptibility with the spin lattice relaxation rate from NMR, as predicted 
by the spin-nematic theory, is found in both electron- and hole-doped BaFe 2 As 2 , but not in FeSe. The 
intricate relationship of the nematic susceptibility to spin and orbital degrees of freedom is discussed. 
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1. Introduction 

Electronic nematic phases have recently attracted considerable attention, especially in connection to 
high-temperature superconductivity [1, 2]. For these phases, the term “nematic”, which originally refers to 
a liquid-crystal phase, in which rotational symmetry is broken while translational symmetry is preserved, 
has been borrowed from its original context to describe a symmetry breaking due to electronic effects [3]. 
In the iron-based systems, the stripe-type antiferromagnetic phase[4-7], which occurs in close proximity to 
the superconducting phase, reduces the C 4 rotational symmetry of the high-temperature tetragonal state to 
C 2 and is therefore referred to as Ising-nematic. It is necessarily accompanied by an orthorhombic lattice 
distortion[8]. Simultaneous magnetic and structural phase transitions occur, for example, in underdoped 
Bai_a;Ka;Fe 2 As 2 [9], suggesting that the small lattice distortion (< 4 x 10“^) is simply induced by magne¬ 
toelastic coupling. However, it was observed early on that the (nearly [10]) concomitant magneto-structural 
phase transition of the parent compound BaFe 2 As 2 splits into two well-defined transitions upon Co substi¬ 
tution [11, 12]. Neutron diffraction studies demonstrated that the structural transition temperature Tg is 
several K higher than the magnetic transition temperature Tjv [13, 14]. Hence, there is a small region of an 
orthorhombic—i.e., “nematic”—but paramagnetic phase in Ba(Fei_a;Coa;) 2 As 2 , which is, e.g., also observed 
in other transition-metal doped 122-systems [15], in F-doped 1111-compounds [16] and in Co-doped NaFeAs 
[17]. Intuitively, Tg > suggests that the structural instability is the primary one, even though a scenario 
in which the structural transition is nevertheless a consequence of magnetic interactions has been proposed 
soon after[18, 19]. The observation certainly puts the simple picture of the orthorhombic distortion as a 
mere consequence of stripe-type antiferromagnetism into question and has sparked great interest in the 
nematic phase of the iron-based materials, including an intensive debate about its microscopic origin [2]. 
Of particular interest in this debate is FeSe, which undergoes a similar structural distortion as the other 
compounds, but does not order magnetically [20]. 

Here we review and discuss our recent measurements of the elastic shear modulus, which is shown to 
be a particularly sensitive probe of the incipient structural distortion, using a novel three-point bending 
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setup in a capacitance dilatometer[21, 22]. We use the notation of the two-iron (or the tetragonal) unit 
cell, in which the structural distortion occurs in the B 2 g channel and is related to the elastic shear modulus 
Cge- The contribution is organized as follows. In the remaining part of this Introduction, we review 
previous experimental and theoretical work on nematicity, in particular as related to the softening of the 
shear modulus above the structural transition. In the following, we detail how Cee can be linked to the 
nematic susceptibility using Landau theory (Section 2). In Section 3 we present briefly the three-point- 
bending technique and in Section 4 we show how the elastic data is influenced by domain formation in 
the orthorhombic state using dynamical three-point bending measurements on BaFe 2 As 2 . In Section 5 we 
present comprehensive shear-modulus data of Ba(Fei_a;Co 2 ;) 2 As 2 and (Bai_a;K2,)Fe2As2[21] and of FeSe[23]. 
Analysis of these data shows that, while the Ba(Fei_a;Coa;) 2 As 2 system can be described in a quantum critical 
scenario, (Bai_a;Ka;)Fe 2 As 2 is characterized by a first-order transition between different ground states on 
increasing K content. The inferred nematic susceptibility of FeSe is shown to be remarkably similar to the 
one of underdoped BaFe 2 As 2 . In Section 6 the nematic susceptibility, as determined by elastoresistivity 
measurements and electronic Raman scattering, is compared to our elastic data. In Section 7 we test the 
scaling of the nematic susceptibility with magnetic fluctuations, derived from NMR data, within the spin- 
nematic scenario, in which the structural transition is a direct consequence of strong magnetic fluctuations. 
This scaling is found to be well satisfied in both Ba(Fei_2;Coa;)2As2[24] and (Bai_a;K 2 ,)Fe 2 As 2 , but fails for 
FeSe. Finally, we present a summary and outlook in Section 8. 

1.1. Electronic in-planc anisotropy 

The anisotropic properties of the magnetic/nematic phase have been studied by various different exper¬ 
imental techniques. An early observation of electronic in-plane anisotropy in the orthorhombic phase was a 
large a—b anisotropy of the resistivity in Ba(Fei_a;Coa;) 2 As 2 (nearly a factor of 2) [25]. However, it was sub¬ 
sequently shown experimentally that this resistivity anisotropy strongly depends on the degree of disorder 
of the samples and the type of substitution [26-31]. Various theoretical works also place emphasis on the 
role of disorder in explaining these observations[30, 32, 33]. Strongly anisotropic defects in the orthorhombic 
state, which are a candidate to produce this resistivity anisotropy, were, indeed, observed using scanning 
tunneling microscopy [34-36] and recent measurements of the anisotropy of the Hall effect also point to a 
dominating role of the carrier mobility anisotropy in creating the resitivity anisotropy[37]. However, optical 
conductivity studies show the importance of orbital anisotropy in addition to the anisotropy of scattering 
rates [27, 38-40]. Recent work stresses, in particular, the importance of the anisotropy of the Drude weight, 
rather than scattering rate, for the resistivity anisotropy[41], as was also indicated by earlier disorder- 
dependent elastoresistivity measurements [42]. Similarly, early angular resolved photo emission spectroscopy 
(ARPES) measurements found a significant shift of the dxz (dyz) orbitals to lower (higher) energies below 
Ts [43]. Further, the thermopower was shown to have an even larger anisotropy than the resistivity in the 
orthorhombic state, arising from an interplay of anisotropic scattering and orbital polarization [44]. Nuclear 
magnetic resonance (NMR) [45] and inelastic neutron scattering studies [46, 47] show significant anisotropy 
of the spin dynamics even for Tg > T > T^. Inelastic neutron scattering as a probe of nematicity is described 
in another contribution to this issue [48]. 

Altogether, the observed large electronic anisotropy in the presence of a rather small lattice distortion 
of (5 = (a — b)/{a + b) < 0.4% supports the idea that an electronic order parameter, which has to be 
“nematic” by symmetry, drives the structural transition. To show this, the resistivity anisotropy induced by 
an externally imposed lattice distortion (applied to a sample via a piezo stack) was measured in Ref. [49]. It 
was found that, indeed, the susceptibility of an electronic nematic order parameter diverges on approaching 
Tg from above, in agreement with the assumption that it drives the structural transition. Similarly, a study 
of the electronic relaxation dynamics using femtosecond-resolved polarimetry suggests that nematicity is an 
independent electronic degree of freedom [50]. 

In order to measure the in-plane anisotropy of various physical quantities, the crystals need to be de- 
twinned, which can be accomplished by the application of uniaxial stress along the tetragonal [110] direction, 
crpio] [25, 51] or, in some cases, by a high magnetic field [52, 53]. An earlier review on detwinning and elec¬ 
tronic in-plane anisotropy in the iron-based superconductors is given in Ref. [51]. However, the application of 
^■[iio] s-lso significantly smears out the structural transition and can induce a marked anisotropy of electronic 
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properties even above Tg [25, 54-56]. This strong sensitivity of the system to cr[iio] already demonstrates 
a large nematic susceptibility and can be used to determine it quantitatively. The nematic susceptibility 
has been evaluated by various techniques, including the above described strain-dependent resistivity[49], the 
stress-dependent optical reflectivity[39], the elastic shear modulus[21] and the Raman response [57]. 

1.2. Two theoretical scenarios - spins vs orbitals 

As to which electronic degrees of freedom, spin or orbital, underlie this electronic nematic order param¬ 
eter, two main scenarios have been discussed. The first one[58-60] places emphasis on orbital degrees of 
freedom, in particular the iron dxz and dyz orbitals, which are degenerate in the high-symmetry tetragonal 
phase. The structural transition, in this scenario, occurs when these orbitals order and become inequivalent 
in energy. The nematic order parameter ip is then given by the difference in orbital occupation. Orbital 
order may trigger a, secondary, magnetic transition [58]. This model naturally explains why the structural 
transition occurs at a higher temperature than the magnetic transition, yet a magnetic transition does not 
necessarily follow the structural one within this orbital scenario. 

On the other hand, magnetism is essential in the second model, where spin fluctuations are considered as 
the driving force for the structural transition. In this spin-nematic scenario[3, 18, 19], the primary instability 
is that of the stripe-type magnetic phase, with the ordering vector eitherQi = (ti'jO) or Q 2 = (0,7r), in 
orthorhombic notation. This state has an additional degree of freedom with respect to, e.g., checkerboard- 
type antiferromagnetism, namely the orientation of the ferromagnetic stripes (i.e., whether the ordering 
vector is Qi or < 52 )) and this is determined by the “spin-nematic” order parameter tp [19]. The calculation of 
Ref. [19] has shown that spin fluctuations can induce a finite ip, and thus a finite orthorhombic distortion, 
at a higher temperature than the magnetic ordering temperature, explaining the observed sequence of phase 
transitions. The debate about whether orbital or magnetic degrees of freedom are “in the driver’s seat” of 
the structural transition has been reviewed in Ref. [2]. 

1.3. Extreme cases: orthorhombic distortion without magnetic order in FeSe and tetragonal magnetic phase 
in hole-doped BaFe 2 As 2 

In most iron-based materials, the onsets of magnetic ordering and the structural distortion occur very 
close to each other in the phase diagram, which indicates that they are strongly coupled and hampers the 
determination of the “driver”[2]. There are, however, two exceptions to this which have recently attracted 
a lot of attention. The first exception is the 11-type iron-based superconductor FeSe, which undergoes a 
tetragonal-to-orthorhombic structural phase transition at Tg ~ 90 K, similar to that found in the underdoped 
122 materials and with similar magnitude of the orthorhombic distortion [20, 23, 61]. Yet, whereas this 
transition always occurs in proximity to stripe-type antiferromagnetic order in the 122 systems, no static 
magnetism was found in FeSe at ambient pressure [20, 62]. However, an enhanced spin-lattice relaxation 
rate in NMR at low temperatures indicates the presence of strong spin fluctuations [63]. In contrast to 
the 122-type systems, phase-pure superconducting FeSe cannot be obtained by self-flux growth techniques. 
Recently however, high-quality single crystals were obtained using low-temperature vapor transport[64] or 
KCI/AICI 3 flux[65, 66 ], enabling a multitude of experimental studies including those of Young’s modulus, 
NMR, elastoresistivity, ARPES and inelastic neutron scattering [23, 67-73]. 

FeSe has a rich temperature-pressure phase diagram [61, 74-76], which has recently been obtained in 
new comprehensiveness using resistivity measurements [77, 78]. The structural transition is found to be 
suppressed by hydrostatic pressure at Y 2 GPa. On the other hand, the low-temperature spin fluctuations 
are enhanced by pressure [63] and probably static magnetic order sets in above ~ 1 GPa with a subsequent 
increase of Tjv [62, 75, 77]. T^, has a non-monotonic pressure dependence, increasing initially and reaching a 
local maximum at ~ 0.8 GPa followed by a slight decrease [75-77]. At higher pressures, increases again 
and the onset of Tc reaches surprising 37 K at ~ 7 GPa [61, 74, 79, 80]. Recent inelastic neutron powder 
[72] and single crystal [73] diffraction results suggest that magnetic fluctuations at ambient pressure occur 
around the same (tt, 0) stripe-type wave vector as in the other iron-based compounds. These puzzling results 
have resently attracted considerable attention and have resulted in various theoretical scenarios[81-85]. 

The second exception where magnetic order and orthorhombic distortion do not closely follow each other 
is the C' 4 -symmetric reentrant magnetic phase in Na-doped BaFe 2 As 2 [ 86 ]. Within a certain range of Na 
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content, the usual orthorhombic distortion first develops with the stripe-type magnetic order below Tn, 
but then suddenly disappears within experimental resolution upon entering this phase. This observation 
was taken as evidence that magnetic degrees of freedom drive the structural phase transition (the spin- 
nematic scenario) because the existence of such a phase can hardly be reconciled with orbital order being a 
prerequisite for magnetism[86]. Yet, the detailed magnetic structure is still under intense study. Polarized 
neutron scattering indicates a spin-reorientation from in-plane to c-axis oriented and suggests that the 
magnetic structure might still be orthorhombic, while a truly tetragonal “2 — Q” magnetic structure could 
not be ruled out [87]. Group theoretical analysis suggests that the question of who is in the driver’s seat may 
be solved by determining the pattern of magnetic and orbital order within this phase [88]. Further theoretical 
works [89, 90] show that interplay between the two stripe-type magnetic ordering vectors Qi = (0,7r) and 
Q 2 = (t!", 0) can lead to the observed phase diagram. This putatively tetragonal magnetic phase was also 
reported in Na-doped SrFe2As2[91] and in single crystals of K-doped BaFe2As2[92], where the interplay with 
superconductivity was also studied. Measurements of the orthorhombic order parameter using capacitance 
dilatometry [92] of these samples yield an upper limit for the orthorhombic distortion of ~ 0.01 x 10“^ in 
the new magnetic phase of Bai_a;K 2 ,Fe 2 As 2 , i.e., less than 1% of the value of S in the stripe-type magnetic 
phase. 

1.4- The elastic shear modulus - soft mode of the structural transition 

It is clear that once a material undergoes the structural transition, all properties (e.g., lattice constants, 
orbital occupation or spin fluctuations) aquire in-plane anisotropy, which makes it difficult to distinguish 
between the two scenarios. As an alternative, the study of the susceptibility of the various quantities above 
the structural transition has been suggested as a possible viable approach to determine the driving force of 
the structural transition [2, 24, 49]. 

In proximity to the structural phase transition of the iron-based materials (which can be classified as 
pseudoproper ferroelastic [49, 93]) the elastic shear modulus Gee is expected to grow soft. Note that the 
inverse of an elastic modulus can be considered a structural susceptibility. The shear modulus Cqq was 
studied extensively in the Ba(Fei_a;COa;) 2 As 2 system using ultrasound [19, 60, 94-96]. These measurements 
confirm Cqq as the soft mode of the structural transition, as expected from the symmetry of the lattice 
distortion. An ultrasound study on single-crystalline FeSe [97] shows the same soft mode. Further, signatures 
of a behavior which was termed “structural quantum criticality” were found around optimal doping in the 
Ba(Fei_a;Coa;) 2 As 2 system. Namely, the temperature and doping dependence of the structural susceptibility 
See, = Cqq was found to resemble closely the magnetic susceptibility in proximity to a magnetic quantum 
critical point [95]. The elastic softening persists over a large part of the superconducting dome on the 
overdoped side which makes the associated fluctuations a possible candidate for the pairing glue in the iron- 
based systems [21, 59]. A pronounced hardening of Gge below the superconducting transition temperature 
Tc in optimally doped Ba(Fei_a;Co 2 ;) 2 As 2 was shown to reflect the competition between magnetism and 
superconductivity because the structural distortion and magnetism are strongly coupled in this system[19]. 

2. Landau theory with bilinear strain-order parameter coupling 

In the Landau theory of a second-order phase transition, a structural distortion can be induced by a 
bilinear coupling in the free energy between elastic strain e and the primary (possibly electronic) order 
parameter ip. In the case of the iron-based materials, the structural distortion is given by 5 = {a—b)/{a + b), 
or £g = -|- the relevant component of the elastic strain tensor, and (p is the electronic nematic order 

parameter. Note that the bilinear coupling between Sg and p is, indeed, only allowed if p is “nematic”, i.e., 
it breaks the four-fold rotational symmetry of the high-temperature phase. However, this phenomenological 
approach is independent of the microscopic origin of p (spin or orbital). Due to the bilinear coupling, Sq 
is a measure of p, i.e., the two quantities are proportional to each other [19]. The stress-strain relation 
o’e = Gggffg in the tetragonal or orthorhombic system shows that the soft elastic mode is Ggg, i.e., a second- 
order structural phase transition occurs when Ggg —> 0. 
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Figure 1: (a) Temperature dependence of the nematic susceptibility Xv? = — ^o) (solid line) and the nematic susceptibility 

renormalized by bilinear coupling to the lattice X‘^ — ~ (dashed line) in a mean-field model, (b) Temperature 

dependence of the soft elastic mode Cee = (^ ~ Tp^)/(T — To) of the structural transition induced by this bilinear coupling 
between the strain component eg and the nematic order parameter eq. 1. The effect of the bilinear coupling is to increase the 
transition temperature from To to Tp^. (c) shows the temperature dependence of Cqq for a range of parameters A^/aC66,0 = 
— Tq. The slope of Cqq just above Tg is determined by the bilinear coupling strength A. 


The relation between the nematic susceptibility and Cqq can be obtained by considering the free energy 
density [8, 19, 49] 

P = Po + 2^66,o£e — ^£6^^ + 2 (Xv) ^( 1 ) 

Here, —X£e(p is the bilinear coupling term (with the coupling constant A), is the bare elastic energy 

{i.e. without the coupling A) and Cee.o is the bare elastic constant, which, by assumption, has no strong 
temperature dependence. The last two terms represent a Landau expansion in the nematic order parameter 
ip, with Xip the bare nematic susceptibility and B > 0 the usual quartic coefficient of the Landau expansion. 
A second bilinearly coupled (nematic) order parameter can be included analogously to describe, e.g., the 
interplay between elastic, spin-nematic and orbital degrees of freedom, as in Refs. [98, 99]. Here, however, 
we restrict ourselves to only one nematic order parameter. In general, the effective, renormalized elastic 
constant Cee is given by [100, 101] 


(fF _ d'^F / d'^F y 

del \de&dip) ydip"^) 


( 2 ) 


because the condition that F is minimal with respect to both eg and ip couples the two order parameters. 
This results in [8, 49] 


Cee — C'ee.o ~ 


A 2 


iXv) V 35(^2 


( 3 ) 


which reduces to Cgg = C'gg.o ~ ^^Xv above the ordering temperature when p = 0. Note that Xv is also 
renormalized by the coupling A and the effective nematic susceptibility Xv is given by {x^p) ^ = {Xip) ^ ~ 
A^/Cgg.o- Rewriting Cgg in terms of Xp^ yields the expression of Ref. [19] 


1 1 A^ ^ 

f-i A f-i‘1 Xp 
>-^66 '-^ 66,0 '-^ 66,0 


( 4 ) 


Note that Eq. 4 is not limited to a Landau expansion in p and is actually valid more generally, as long as 
bilinear coupling to a harmonic lattice is considered [19]. Equations 3 and 4 show that we can access the 
nematic susceptibility Xp s-nd also Xp hy measuring Cgg. 

Before analysing real data, we illustrate the expected behavior of the shear modulus in the above Landau 
theory. We assume a mean-field Curie-Weiss-type divergence of Xp approaching Tq, which would be the 
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transition temperature in the absence of coupling to the elastic strain, 


a(r-To)' ^ ’ 

This leads to the temperature dependence of the soft elastic mode in a mean-field case [101], 

Cee = Cee.o j for T > (6) 

f 2(Tf^-T) \ f.. 

o-ee — 0/66.0 ^ ^j^cw _j'^_2T J ^ ^ ’ I'l 

shown in Fig. 1, with = Tq + ^ the new transition temperature, increased with respect to Tq 

by the coupling to the strain. An energy scale characteristic of the coupling is given by — Tq, and 

determines the curvature of Cqq{T) in Fig. 1 b, c. Note that the renormalized = l/a{T — T^^) naturally 
diverges at the new transition temperature. 

This above model, however, is limited to a description of a second-order nematic/structural transition 
and neglects the subsequent (or concomitant) magnetic transition occurring in most iron-based systems. To 
illustrate the expected temperature dependence of Cqq for BaFe 2 As 2 , the Landau model can be expanded 
to describe this split magneto-structural transition phenomenologically. This is achieved by the free energy 
density 


F = Fo+^{T- To)(p2 + ^ _ Tn,o)M^ + (8) 

which is equivalent to equation 1 concerning the elastic and nematic contributions and adds a magnetic 
order parameter M. Additionally, the contribution of conjugated uniaxial stress —Sea is also included to 
model the behavior under finite stress. For parameters values a = 1, Tq = 1.025, B = 1, C66,o = 1, A = 0.2, 
M = 1, Tjv^o = 1, f = 1 and /i = 0.2 the model reproduces the second order structural and a first-order 
magnetic transition slightly below the structural one. Note that the coupling between M and ip causes the 
magnetic transition to be first order despite a positive fourth-order coefficient v. A solution which minimizes 
F is calculated numerically for varying values of cr and is shown together with the resulting elastic modulus, 
C66 (still given by equation 3) in Fig. 3b. Note that the temperature dependence of Cee is unchanged with 
respect to the solution of eq. 1 in the high-temperature region where M = 0. 

In principle, the structural transition may arise either from a divergence of Xvi s-s in an electronically- 
driven transition, or from a vanishing of C 66 , 0 ) s-s in a bare lattice instability. It has been proposed that the 
temperature dependence of Cqs close to the phase transition can distinguish the two cases [8], namely C66 
should vanish linearly in the case of a bare elastic instability or with a curvature as in Fig. 1 in the case of 
an electronically driven transition, which is assumed here. 

We note that the above Landau formalism is a mean-field treatment, in which fluctuations of the order- 
parameter are ignored. In fact, this often turns out to be a good approximation for second-order structural 
phase transitions, because the long-range elastic interactions strongly suppress fluctuations [102-104]. 

3. Technique: Three-point bending in a capacitance dilatometer 

Three-point bending is a long-standing and widely used mechanical test to study elastic properties of 
diverse materials. The technique is particularly appealing by its simplicity. A platelet- or beam-shaped 
sample is supported along two lines, while a force is applied at a third, middle, line and the induced 
deflection is measured. Small forces result in a sizable deflection so that the elastic modulus can be measured 
comparatively easily. In the limit of thin samples, the measured elastic modulus is proportional the material’s 
Young modulus, which also determines the sample stiffness in a uniaxial tension/compression experiment 
[105], even though samples are subject to non-uniform stress in three-point bending[106]. We have developed 
a three-point bending setup in a capacitance dilatometer[107] to measure the Young modulus of iron-based 
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Figure 2: (a) Schematic representation of the capacitance dilatometer with a sample inserted for three-point bending. The 
sample (see panel (b), of dimensions Zxtux/i^SxlxO.l mm®) is supported by three wires and pressed against one movable 
plate of a plate-type capacitor using a screw with a force of ^ 0.2 N. The movable capacitor plate is suspended via a set of 
circular parallel springs and the change of the capacitor gap d (indicated by arrows) on changing temperature is measured, 
(c) shows the regular setup with a sample inserted for thermal expansion measurement, (d) Temperature evolution of the 
capacitor gap d{T), expected in a c-axis thermal-expansion (TE) experiment of a O.f mm thick BaFe 2 As 2 sample[108] (black 
dashed line), and measured in three-point bending (3PB) of such a sample oriented along [100] (blue line) and [110] (red line). 
An exceedingly large effect is observed in the latter case, a result of the strong temperature dependence of the shear modulus 
Cee of BaFe 2 As 2 . 


materials with very high resolution [21]. The advantage of this technique over the ultrasound measurements 
used more traditionally to determine elastic constants, is that the sample requirements are not as stringent. 
In particular, ultrasound velocity measurements require samples of considerable size and quality. On the 
other hand, the typically quite thin platelet-like single crystals available for the iron-based based materials 
are perfectly suited for the three-point bending technique. This has allowed us to investigate a larger variety 
of materials than would have been possible using ultrasound. 

In a typical capacitance dilatometer, the sample is held in place by a small uniaxial force from the two 
leaf springs of the parallelogram arrangement holding the moveable capacitor plate (see Fig. 2c). Hence, 
sample-length changes (due to, e.g., thermal expansion) lead to a change of the capacitor gap d. A very high 
length resolution of 0.1 — 0.01 A can be achieved by measuring the resulting change of capacitance [107]. By 
placing a sample in three-point bending configuration (Fig. 2b) into the dilatometer, as shown in Fig. 2a, 
one no longer measures the thermal expansion of the sample, but rather its elastic bending modulus. This 
is because the small force from the dilatometer causes the sample to bend and the effective height of the 
arrrangement is determined by the sample’s flection, hence, its elastic modulus. Notably, when the elastic 
modulus of such a sample becomes soft (i.e., decreases), the sample bends more strongly so that its effective 
height along the axis of the dilatometer decreases. Fig 2 (d) shows examples of such measurements. The red 
and the blue curves show three-point bending experiments with samples oriented with their tetragonal [110] 
and [100] directions, respectively, perpendicular to the supports. An exceedingly large effect is observed in 
three-point bending along [110]. This is because the Young’s modulus along [110], 
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(9) 


is dominated by the elastic shear modulus Cge as long as Cqq is smaller than the other , and Cqq decreases 
strongly on cooling towards Ts[95]. In contrast, the black dashed line shows the much smaller temperature 
dependent change of the capacitance gap Ad expected from just the thermal expansion of a 100 /im thick 
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Figure 3: (a) Real part of the Young modulus BaFe 2 As 2 obtained in dynamic three-point bending experiments in which 

the applied force F = Fs + Fj^) ex.p(iujt) varies at a frequency of a;/27r = 1 Hz (see inset). Ypig] ^ ^dI'^d exp(f(5) is calculated 
from the measured displacement u = us -\-U£) exp{i{ut — (5)) under the applied force (see inset). Data are taken at different Fs 
and constant ratio Fjj/Fs- The maximal uniaxial stress at the sample surface arising from Fs, crmax^S = 2 wh^ [106] is given 
as a characteristic parameter. The two measurements with the lowest force were conducted in “constant amplitude” mode, 
in which Fs and Fd are smaller when the sample is soft and increase when the sample becomes harder, thus increasing the 
resolution, (b) Soft elastic mode in a mean-field picture for closely spaced structural and magnetic phase transitions, obtained 
via numerically minimizing eq. 8 (see text for parameters) and then calculating Cqq. Good agreement with the experimental 
data is obtained for the high-stress measurement, while the measurements at lower stress and below Ts are dominated by the 
effect of structural twins. Insets show the order parameters </? and |M| of the numerical solution. 


sample of BaFe2As2[108]. A detailed description of the quantitative analysis of the data is given in [22]. 
The technique is particularly well suited to access large changes of the Young modulus of a sample occurring 
over a broad temperature range, as well as any sharp anomalies. As we will show, this technique even has 
the resolution to detect the often very small anomalies at the onset of superconductivity. 

In addition to the static three-point bending measurements in the capacitance dilatometer, ultralow- 
frequency dynamic three-point bending measurements of BaFe 2 As 2 were performed using a Diamond DMA 
(dynamical mechanical analyzer) from PerkinElmer [109]. Here, samples are subjected to a force F = 
Fs + Fd exp(iu}t) with a dynamic component at the very low frequency of 1 Hz. From the measured 
displacement u = us + ud — 5)), the dynamic stiffness ks = Fd/ujj expiiS) oc Y[iio] is calculated. 

kg may be a complex number in the presence of dissipation, and its real part corresponds to the results of 
the static measurements. 


4. Dynamical three-point bending measurements of BaFe 2 As 2 

Before addressing the detailed static measurements of the Young’s modulus in a capacitance dilatometer, 
we show in this section results of dynamic three-point bending experiments on BaFe 2 As 2 , which elucidate 
the role of structural twins in the ordered phase[110]. For these measurements the sample is subjected to the 
force F = Fs + Fjj exp(iwt) with a dynamic component Fjj ^ (0.7 — 0.8)Fs at a frequency of lij/2tt = 1 Hz. 
Fig. 3a shows the real part of the Young modulus (equivalent to the result of static measurements) 

for different Fs but constant ratio F^/Fs- 

A smooth softening of Fjno] upon cooling towards the structural transition is observed (Fig. 3a), in 
agreement with the soft-mode behavior shown in Fig. 1. However, at the lowest stress, a kink and a 
small, sharp, minimum may be ascribed to the structural and magnetic transition, respectively, and Y[iio] 
is essentially temperature independent below T/Vj in strong contrast to the behavior shown in Fig. 1. With 
increasing stress, the data deviate from the low-stress curve below some T > Tg and Fjno] does not soften as 
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much in total. The structural transition is smeared out, while the magnetic transition is affected much less. 
Notably, a discontinuous hardening of the Young modulus on cooling through T/v emerges with increasing 
stress. It is reminiscent of the phonon modes measured using inelastic neutron scattering [111]. At the 
highest stress, Tfug] indeed recovers its initial high-temperature value at ~ 0.9Tn. Note that the deviation 
of this curve from the others at higher temperatures is an experimental artefact. 

The experimental data at high stress are well reproduced by the model in eq. 8 of a split magneto- 
structural transition (Fig. 3b). However, the experimental data obtained at lower stress are much too 
small in comparison to this model for temperatures T <Tg. These anomalously low values of the Young’s 
modulus below Tg likely arise from the ’superelastic’ behaviour of the twinned samples, which is due to 
domain wall motion and typically found in ferroelastic materials[105, 112, 113]. The effect can be understood 
by considering two kinds of structural domains, '-|-^-type domains which are elongated along [ 110 ] and ' 
type type domains which are shortened along this direction. Application of a small compressive stress 
(T[iio] increases the fraction of '—'-type domains at the expense of the '-|-'-type domains, so that the total 
sample length decreases through domain wall motion. A relatively small applied force can thus induce 
significant length changes, which means that the elastic modulus appears very small at low frequencies. 
In the three-point bending experiment, samples are essentially subjected to stress along the [110] direction, 
which suppresses domain formation and leads to a strong variation of the measured Fjno] with applied stress. 
The true monodomain elastic properties of the system can only be obtained in these dynamic three-point 
bending experiments if the applied stress is strong enough to force the system into a single domain state. 


5 . Young’s modulus of Ba(Fe,Co)2As2, (Ba,K)Fe2As2 and FeSe 

Figure 4 shows the Young’s modulus Y[iio] of Co- and K-doped BaFe 2 As 2 and of FeSe, as measured 
by the three-point bending technique in a capacitance dilatometer described in section 3 [21, 23]. Fig. 4c 
shows that the temperature dependence of Ypig] obtained by three-point bending is very similar to Cge 
(see eq. 9) as determined by ultrasound measurements for Ba(Fei_a;Coa;) 2 As 2 [95], the system for which 
such data is available. Notably, all essential features of Cgg are also observed in the Young’s modulus 
data, confirming the capacity of our experiment to determine the temperature dependence of the shear 
modulus. In BaFe 2 As 2 , strong softening of Yj^^gj upon approaching the transition at Tg from above is 
observed. The softening gradually disappears with both Co- and K-doping. Interestingly, the Young modulus 
Y[iio] of FeSe, which undergoes a similar structural transition as the iron-arsenides but no magnetic phase 
transition, shows the same pronounced softening on cooling towards « 90 K as underdoped BaFe 2 As 2 
(Fig. 4(d)). For moderately overdoped samples, a marked softening on cooling is still observed, while 
Y[]^ 20 ] hardens considerably below Tc- The Young’s modulus of a non-superconducting, strongly overdoped 
Ba(Feg, 67 Coo, 33 ) 2 As 2 sample is found to harden by a few % on cooling from room temperature, which is the 
typical behavior induced by phonon hardening [114]. The fiat temperature dependence in the orthorhombic 
state below Tg is ascribed to the multidomain effect outlined in the previous section and, in the following, we 
will therefore focus on data at T > Tg. We note that we find no evidence of a higher-temperature {T > Tg) 
nematic transition as proposed by Kasahara et ah, [115] in any of our high-resolution shear modulus data. 
This is in agreement with previous high-resolution thermal-expansion measurements on both Co- and P- 
substituted BaFe 2 As 2 [ 108 , 116] and heat-capacity measurements in BaFe2(Asi_2,Pa;)2[117], which also do 
not find evidence of an additional phase transition. 


5.1. Doping and temperature dependence of the nematic susceptibility 

Eq. 3 directly links Cqq to the nematic susceptibility Xv under the assumption that the structural 
transition is induced by bilinear coupling to a nematic order parameter, eq. 1. In real systems the bare 
elastic constant < 766,0 is slightly temperature dependent because of phonon anharmonicity and typically 
hardens by a few % upon cooling between room temperature and zero temperature [114]. In order to 
remove this background contribution and to connect with the Young modulus, we make the approximation 


Y[iio] ^ <^66 
Yg (766,0 


( 10 ) 
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Figure 4: Young’s modulus along the tetragonal [110] direction Yjiio] of (a) Ba(Fei_a;Coj;) 2 As 2 and (b) Bai_a;Kj;Fe 2 As 2 
for a wide doping range, measured using the static three-point bending in a capacitance dilatometer described in section 3. 
The data for T < Tg are dominated by the effect of structural twins, similar to the low-stress measurements of Fig. 3. (c) 
Young’s modulus (lines) and shear modulus Cqq (symbols, from Ref. [95]) for Co-doped BaFe 2 As 2 . Both quantities 

are normalized at room temperature. The indicated Co content is taken from Ref. [95] and samples with similar transition 
temperature are compared. All essential features, namely the strong softening of Yjno] cooling towards Ts, and the marked 
hardening below Tc for moderately overdoped samples, are reproduced in the Young’s modulus data. Those curves show, 
however, a slightly stronger curvature than Cee, presumably because of the contributions from the other elastic constants (see 
eq. 9). (d) Young’s modulus Yj^o] single crystalline FeSe, which shows very similar softening on cooling towards T^. 
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Figure 5: Nematic susceptibility x<^ units of 066,0/-^^ of {^) Ba(Fei_a;Cox) 2 As 2 and (b) Bai_a;Ka;Fe 2 As 2 obtained from 
the data in Fig. 4. For a second order structural phase transition should reach 1 at Ts and it is unclear why the 

experimental data do not reach this value, (c), (d) show the inverse revealing a Curie-Weiss-like temperature 

dependence of the nematic susceptibility for all samples except Bai_xKccFe 2 As 2 with > 30% K content. The dashed straight 
line in (c) is a guide to the eye. (e) and (f) show similar data for FeSe, which resemble closely underdoped BaFe 2 As 2 . 
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Figure 6: (a) Nematic susceptibility as a color-coded 

map in the composition-temperature phase diagram of 
Ba(Fei_j;Coa;) 2 As 2 and Bai_j;Kj;Fe 2 As 2 . Green trian¬ 
gles show the inflection point of (T) as a lower limit for 

the validity of Curie-Weiss law \'^X<p/^66,0 = ^ 

(see text), (b) Weiss temperature To and temperature at 
which 6*66 extrapolates to zero, Tp^. When the data 
deviate from the simple Curie-Weiss law, as for samples 
with > 30 % K content, this deviation is ascribed to a 
temperature dependent To and To (extrapolated) at var¬ 
ious temperatures is reported. The dashed area indicates 
the doping range, in which an abrupt change suggests a 
first order transition between different ground states, (c) 
Curie constant A^/aC66,0 = — Tq. The dashed line 

shows the extrapolation for the region where the data do 
not follow the Curie-Weiss law, used to obtain To(T) in 
(b). 


Here, Yq is the non-critical contribution to for which we use the data from a 33% Co substituted 

sample and further assume that it is independent of doping[108]. Additional details can be found in [21, 22]. 
Assuming that the structural transition of FeSe is also induced by bilinear coupling between a nematic order 
parameter and the lattice, (p is inferred from the Young’s modulus data of FeSe in the same way, assuming 
similar Yo[23]. The obtained bare nematic susceptibility is plotted in Fig. 5 (a), (b), (e) in units of 
A^/Cee.o- Note that the nematic susceptibility renormalized by the coupling to the lattice can be inferred 
easily using (A^x^^/Cee.o) ^ = (A^Xip/Cee.o) ^ ~ 1- X< 4 > increases strongly with decreasing temperature for 
all but the most strongly overdoped BaFe 2 As 2 samples. Interestingly, FeSe is found to resemble closely 
underdoped BaFe 2 As 2 . 

The inverse {x>^) plotted in Fig. 5 (c), (d), (f), evidences a mean-field-type temperature dependence 
Xv = — Tq)) for Ba(Fei_a;Co 2 ) 2 As 2 up to x = 0.09, in agreement with Ref. [95], for Bai_a;Ka;Fe 2 As 2 

up to X = 0.24 and for FeSe. However, and somewhat surprisingly, Ffno] does not follow a Curie-Weiss law 
in the higher doped Bai_a;K 2 ,Fe 2 As 2 samples in which the structural transition is suppressed. Instead, the 
temperature dependence is found to be less “critical” and Yjhq] displays a clear inflection point. The detailed 
doping and temperature dependence of Xip is presented in Fig. 6. The color map of Xv i’^ Fi§- shows that 
the nematic susceptibility is enhanced fairly symmetrically over most of the superconducting domes of both 
Ba(Fei_a;Coa;) 2 As 2 and Bai_a;K 2 ,Fe 2 As 2 , suggesting its possible role in promoting superconductivity[21]. To 
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describe the temperature dependence of Xv for different doping levels, we consider two parameters. First, 
we define T* as either the infiection point of or its maximum value, whichever is greater. T* is thus a 

lower limit for the validity of a Curie-Weiss law. Second, deviations from the Curie-Weiss law are ascribed 
to a temperature dependence of the parameter Tq, whose value at fixed temperatures is plotted in Fig. 6 b. 
T* reaches near zero temperature around optimal doping, consistent with a quantum critical scenario [95] 
only in the Ba(Fei_a;Coa;) 2 As 2 system. In contrast, T* does not go below ~ 75 K in Bai_a;Ka;Fe 2 As 2 . The 
findings suggest a first-order transition between different ground states preempting a quantum critical point 
in Bai_a;Ka;Fe 2 As 2 , in particular the step-like change of Tg between 24% and 30% K content (Fig. 6 b). 
Interestingly, a new, C 4 -symmetric, magnetic phase was subsequently found to emerge in between these two 
doping levels [92]. 

In Ref. [118], Cqq{T) was calculated including vertex corrections and was found not to follow the mean- 
field-type Curie-Weiss law. The data in Figs.4 and 5 can be well described in this approach, for which 
the change of behaviour of Bai_a;Ka;Fe 2 As 2 between 24% and 30% K content would be a consequence of a 
simple change of parameters as might be due to, e.g., a Lifshitz transition[118]. Finally, it is notable that 
the Curie constant A^/aCge.o = — Tg « 30 — 40 K, which is the characteristic energy of the electron- 

lattice coupling, is found to be nearly doping independent up to 9% Co content and 24% K content, the 
compositions for which it can be evaluated. Note that a value of ~ 50 — 60 K is found for same parameter by 
evaluating Cgg ultrasound data in the Ba(Fei_a;Coa;) 2 As 2 system [95]. This value is comparable but slightly 
larger than our result and refiects the lower curvature of the Cgg with respect to the Ffugj data in Fig. 4c. 

5.2. Young’s modulus around Tc 

Besides the softening due to the nematic transition, the high resolution of our three-point bending 
experiment can also be used to study the behavior at and below the superconducting transition. In Fig. 7, 
we show our Young’s modulus data for Co- and K-doped BaFe 2 As 2 as well as for FeSe around Tc in detail. 
For most samples, a pronounced hardening of Y[iig] below Tc is observed, while for some compositions of 
strongly overdoped BaFe 2 As 2 and FeSe a small step-like softening of F[iic] at Tc is also visible. A hardening 
of the elastic shear modulus Cgg below was reported previously in overdoped Ba(Fei_a;COa;) 2 As 2 [19, 95]. 
As explained below, this effect is unusual in its sign, shape and magnitude for the effect of superconductivity 
on an elastic modulus. 

The usual thermodynamic signature of a second-order phase transition in the Young’s modulus or any 
other elastic modulus is a small sudden softening upon cooling, which is related to the strain/stress depen¬ 
dence of Tc. Such a discontinuity is expected because the Young’s modulus is the inverse of a component 
of the elastic compliance Sa, which, in turn, is a second derivative of the free energy. An Ehrenfest-type 
relation relates the size of this discontinuity at Tc, AYi, to the uniaxial-pressure derivative dTc/dpi, 

( 11 ) 

Here, ACp > 0 is the specific-heat discontinuity at Tc and i stands for the direction. Note that necessarily 
AYi < 0, i.e., Yi shows a step-like decrease on cooling through the transition. Also, the shear modulus Cgg 
alone is not expected to have such discontinuity, because the first derivative of Tc with respect to a shear 
deformation necessarily vanishes [119]. Hence, any discontinuity in the Young’s modulus arises from the 
contribution of the longitudial elastic constants (the ’ 7 ’ in Eq. 9) and AYpigj/Yj^^Pj = AYjiggj/Yj^Ppj. It has 
been pointed out in Ref. [120] that for an exotic superconducting state which mixes nematic fluctuations and 
s- and d-wave superconductivity and breaks tetragonal symmetry, even Cgg can have a discontinuity at Tc. 
Such a state might occur in overdoped Bai_ 2 ,Ka;Fe 2 As 2 [120], however, the Young’s modulus is not an ideal 
probe to search for this effect, since it also contains the contribution from longitudinal elastic constants, as 
mentioned above. 

Fig. 7 presents our Young’s modulus data close to Tc in detail. The hardening of Yfugj below Tc is 
resolved for Bai_a;K 2 ,Fe 2 As 2 up to 82% K content and Ba(Fei_a;Co 2 ;) 2 As 2 up to 9% Co-content (panels 
(a), (b), (j), (k)). Note that the effect is even observed for underdoped samples, when Tc lies within the 
orthorhombic state and the Young’s modulus is actually dominated by the effect of structural twins (see 
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Table 1: Uniaxial in-plane pressure derivative \dTc/dpa\ (sixth column), inferred using the Ehrenfest relation, eq. 11, from the 
discontinuity of the Young modulus at Tc (third column, Fig. 7) and the specific heat (fifth column). A high-temperature 
Young modulus Yi(300 K) = 80 GPa has been assumed in all cases. The obtained value should be the same for the two 
in-plane inplane directions i = [110] and i = [100]. For comparison, the last column shows dTcjdpa inferred from uniaxial 
thermal-expansion and specific heat data. The sign of dTcjdpa can only be obtained from thermal expansion. 



i 

AF, 

U(U) 

ACp/Te 

mJ mol-^K”^ K 

defs. [64, 123, 126] e 

dTa 

dTg 
dp a 

'a K/GPa 

1 Refs. [23, 64, 108] 

Fi(300 K) 

10-3 

F(300 K) 

1 

dp a 

/Gf 
q. 1 

Bao.52Ko.48Fe2As2 

[100] 

-0.62 

1.03 

126 


1.9 


-2.3 

Bao.52Ko.48Fe2As2 

[110] 

-1 

0.92 

126 


2.7 


-2.3 

Bao.4oKo.6oFe2As2 

[110] 

-0.73 

0.99 

96 


2.5 


-2.5 

Bao.ioKo.9oFe2As2 

[100] 

-0.25 

1.07 

50 


1.8 


-2.2 

Bao.ioKo.9oFe2As2 

[110] 

-0.18 

1.05 

126 


1.6 


-2.2 

Ba(Feo,88Coo,i2)2As2 

[110] 

-0.8 

1.02 

9 


8 


« 12 

FeSe 

[110] 

-0.0051 

0.16 

5.6 


3.2 
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section 4). As detailed in Ref. [19], a hardening of the shear modulus Cqq shows that nematic fluctuations, 
which decrease Cee, are weakened in the superconducting state. Its observation up to 82% K content 
indicates the presence of such nematic fluctuations over a large part of the superconducting dome in the 
phase diagram of Bai_a;Ka;Fe 2 As 2 , although the degree of hardening strongly decreases with doping. Since 
this kind of response of V[iio] to superconductivity derives from the contribution of the shear modulus Cge, 
it should not be present in the Young’s modulus along [100], Y^qq], which does not contain any contribution 
from Cqq. Consistenly, panels (c), (d), (e) show that such a hardening is not present in Y[ioo] of overdoped 
Bai_a;Ka;Fe 2 As 2 and very small effects may be due to sample misalignment. 

In FeSe, such a hardening of Yfno] setting in abruptly at Tc is not observed, which can be taken as a sign 
that orthorhombic distortion and superconductivity do not compete in the same way as in doped BaFe 2 As 2 
in this material[23, 64]. However, there is a very slight hardening of (two orders of magnitude smaller 
than for similar Ba(Fei_ 2 ;Coa;) 2 As 2 samples) with an onset well above Tc, which correlates with an anomalous 
uniaxial thermal expansion seen previously [64]. The origin of this effect is unclear at this point and it may 
possibly be a consequence of the presence of very small Fermi energies in the system[69, 70, 121, 122]. For 
some compositions of strongly overdoped Bai_a;K 2 ,Fe 2 As 2 (Fig. 7 (f), (g), (i)), for Ba(Feo. 88 Coo,i 2 ) 2 As 2 
(Fig. 7 (1)), and for FeSe (Fig. 7 (n)), we can also resolve the small step-like softening of Y[iio] and Y[ioo] 
expected from thermodynamics (Eq. 11). Using additional speciflc-heat data [23, 123, 124], the pressure- 
derivative of Tc is calculated via eq. 11 and shown in table 1. The uniaxial pressure derivative dTc/dpi can 
also be calculated by using uniaxial thermal expansion and specific heat via a similar Ehrenfest relation 
dTc/dpi = Vm^Cii/^{Cp/T). As shown in table 1, the pressure derivatives derived from evaluating either 
the Young’s modulus or the thermal-expansion [64, 92, 125] data agree quite well. It is unclear why we do 
not resolve this step-like softening of the Young’s modulus in the Bao.i 8 Ko. 82 Fe 2 As 2 sample nor in Y[ioq] of 
the Bao, 4 Ko. 6 Fe 2 As 2 sample. These, instead, seem to show a step-like hardening that cannot be explained 
with eq. 11. 

6. Different probes of the nematic susceptibility 

Quite generally, the nematic susceptibility can be probed by measuring the sensitivity of the electronic 
anisotropy of various quantities to uniaxial stress or strain [49]. The most detailed data have been obtained 
using strain-dependent resistivity[42, 49, 127, 128]. The Raman response in different symmetry channels 
provides another probe [57] without actually having to apply any stress or strain to the crystal. In this section, 
we briefly review these studies and compare them to the results of our shear-modulus measurements. We note 
that there are also other probes, which were used to obtain the nematic susceptibility, e.g., stress-dependent 
measurements of the optical reflectivity[39, 40]. 
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Figure 7: Enlarged view of the low-temperature Young modulus Yj^io] of (a) overdoped and (b) underdoped Bai_ 2 ;K 3 ;Fe 2 As 2 , 
showing a pronounced hardening below up to 82% K content, (c)-(e) Low-temperature Young’s modulus along [IfO], Yjiio], 
(solid line) and along [100], Y[ioo]i (dashed line) of overdoped Bai—a;K 2 ,Fe 2 As 2 , showing that the hardening occurs only in the 
[110] direction, (f)-(i) Enlarged view of the data close to Tc showing an additional small step-like anomaly AY; at Tc- A linear 
background has been subtracted from the data, (j)-(l) Same as in (a), (b) for Ba(Fei_a;Coa;) 2 As 2 . (m) Enlarged view of the 
low-temperature Young modulus of FeSe. (f) Even further enlarged view of the Young modulus of FeSe around Tc, with a 
linear background subtracted. Vertical arrows mark Tc in all panels. Data in (a)-(e) and (j)-(m) have been shifted vertically 
for better comparison. Lines in (i), (n) indicate an ideal second-order phase transition. 


6.1. Elastoresistivity 

The elastoresistivity is defined as the resistance change induced by sample deformation (strain) and is 
closely related to the piezoresistivity which is the resistance change due to stresses acting on the sample [127]. 
Since the in-plane resistance anisotropy can be taken as a proxy for the nematic order parameter in the Fe- 
based systems, there is an elastoresistivity coefficient, namely mgg, which is directly linked to the nematic 
susceptibility. In Refs. [49, 127, 128], the in-plane resistivity anisotropy N of iron-based materials was 
measured as a function of strain eg, externally applied to the sample via a piezo stack. Making use of the 
bilinear coupling A between nematic order parameter and shear deformation (eq. 1), one obtains 2TOgg = 
diV/deg cx dipjdeQ = XX(p [49, 128]. The proportionality constant between N and ip is related to the details 
of the Fermi surface and electronic scattering and can, in principle, be temperature, as well as doping 
dependent. It was shown that mgg follows a Curie-Weiss law ~ 1/(T —Tg) in BaFe 2 (Aso. 7 Po. 3)2 [128], which 
advocates that this proportionality constant is only very weakly temperature dependent. Further, Ni- and 
Co-doped samples with the same Tg were compared [42]. Ni doping induces larger scattering (lower RRR), 
however, it was found that the elastoresistivity mgg is independent of this kind of disorder for T > Tg [42]. 

Elastoresistivity measurements of the optimally doped compounds BaFe 2 (Aso. 7 Po. 3 ) 2 , Bao.eKo, 4 Fe 2 As 2 , 
Ba(Feo, 93 Nio.o 7 ) 2 As 2 , Ba(Feo. 955 Nio,o 45 ) 2 As 2 , Fe(Teo. 3800 , 4 ) [128] as well as in FeSe[ 68 ] show that the ela¬ 
storesistivity coefficient mgg diverges following an approximate Curie-Weiss law in all of these systems, 
though deviations below ~ 100 K are sometimes observed. We note that the quantity A^y^^/Cgg^o obtained 
from the shear modulus is normalized such that it reaches (ideally) the value 1 at the phase transition when 
Cgg is expected to vanish. In contrast, the elastoresistivity is not normalized and the value of mgg is indeed 
found to be doping dependent. It seems, e.g., to peak around optimal doping in Ba(Fei_ 2 ;COa;) 2 As 2 which 
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was suggested to be due to enhanced fluctuations near a quantum critical point [128]. However, mee still as 
a similar magnitude in all of the studied iron-based systems, although its sign depends on the particular 
system. 

As mentioned previously, the structural transition may, in principle, arise either from a divergence of 
Xip, as in an electronically-driven transition (as assumed here), or from vanishing (766,0) as in a bare lattice 
instability. The quantity obtained by the shear modulus measurements, X'^Xip/C^66,o diverges in both cases. 
Notably, the elastoresistivity rriee oc Xxp is independent of the bare shear modulus (766,o- If the lattice caused 
the transition, then A^x<^/C'66,o would diverge, but would show no strong temperature dependence. The 
experiment by Chu et al.[49] showed that, however, Xx^p also diverges, which was taken as a proof that Lp 
drives the transition and that the lattice distortion is just a consequence of the bilinear coupling to (^[49]. 

6.2. Electronic Raman scattering 

Electronic Raman scattering of Ba(Fei_a;Co2;)2As2[57, 129] finds an enhancement of the Raman response 
in the i? 2 g symmetry channel (which is the symmetry that corresponds to the orthorhombic distortion) 
with respect to the Big symmetry channel. Using Kramer’s Kronig relations, the static nematic charge 

susceptibility Xo ^ extracted from the data. Xo ^ of Ba(Fei_a;Co 2 ,) 2 As 2 was found to increase 

on approaching the structural transition, though it does not diverge at Tg. It, indeed, follows a Curie- 
Weiss law for a wide doping range [57]. Recently, Raman scattering data and Xo ^ also reported for 

FeSe and Bai_a;Ka;As2As2[130]. While Xo ^ of FoSe and underdoped Bai_a;Ka;As 2 As 2 shows the same 
Curie-Weiss-like divergence, deviations from the Curie-Weiss law were observed for close to optimally doped 
Bai_a;Ka;As 2 As 2 samples. 

In Ref. [118], both the shear modulus and the Raman susceptibility have been calculated in the five 
orbital model including vertex corrections. It was noted that the Raman susceptibility is less singular than 
the shear modulus because the photons cannot couple to the acoustic lattice vibrations due to the mismatch 
of their wavelength for the same frequency [118]. In this sense, the charge nematic susceptibility in Raman 
scattering should be similar to the bare nematic susceptibility Xv not Xv) which is renormalized by 
coupling to the lattice[131]. Electronic Raman scattering as a probe of nematicity in iron-based systems is 
discussed in another contribution to this issue [132]. 

There is very good general agreement between the nematic susceptibility from shear-modulus measure- 

2 _ 2 2 _ 2 
ments Xip and the charge nematic susceptibility Xo ^ inferred from Raman scattering. First, Xc/) and Xo ^ 

both follow the familiar Curie-Weiss temperature dependence ~ 1/{T — Tq) [57, 129] in Ba(Fei_a;COa;) 2 As 2 , 
as well as in FeSe[130, 131]. In particular, the Weiss temperature Tq obtained by the two probes agrees very 
well in the Ba(Fei_a;Coa;) 2 As 2 system (Fig. 8, see below). For optimally doped Bai_a;Ka;Fe 2 As 2 , the shear- 
modulus data (see Fig. 5) show an inflection point and strong deviations from a Curie-Weiss dependence of 
the nematic susceptibility. Interestingly, a very similar temperature dependence of the Raman susceptibility 
Xo ^ has also been observed in close to optimally doped Bai_2,Ka;Fe2As2[130]. On the other hand, the 
elastoresistivity coefficient meo of all studied materials, including optimally doped Bai_a;Ka;Fe 2 As 2 , shows 
an approximate Curie-Weiss-like divergence[128]. The origin of this difference is unclear to us. 

Among the above materials, the Ba(Fei_a;Coa,) 2 As 2 system has been studied most intensively. In this 
system, all three Xv) Xo ^ and TO66 show a Curie-Weiss like divergence. To make a quantitative comparison, 
we show in Fig. 8 the values of the Weiss temperature Tq obtained from fitting these data to a Curie-Weiss law 
~ l/(r —To)- Note that in case of the Weo data, the low-temperature region, where meo deviates somewhat 
from a Curie-Weiss law, has been excluded from the flt[128]. We have also included Tq^TiT obtained from 
fitting the spin-lattice relaxation rate divided by temperature l/riT[133, 134], since it also follows a Curie- 
Weiss temperature dependence [133]. 1/TiT is related to the strength of magnetic fluctuations, which the 
spin-nematic scenario links to the nematic susceptibility (these data are discussed in the following section 
7). Also shown in , the temperature at which Cqq extrapolates to zero, i.e., at which the structural 
transition would be expected (see section 2). Why Cqq does not quite reach zero and is lower than Tg 
is still an open question. 

Intriguingly, Tq as determined by electronic Raman scattering, agrees quite well with the elastic data 
over the whole doping range, which supports that the charge nematic susceptibility Xo ^ ol iFe Raman 
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Figure 8: Weiss-temperature To obtained from fitting the nematic susceptibility as determined by Young’s modulus yjiio][21] 
(red circles), electronic Raman scattering[57] (orange squares) and elastoresistivity m66[128] (black open stars) measurements 
in Ba(Fei_a;Coa;) 2 As 2 . Also given is the Weiss temperature of the spin-lattice relaxation rate divided by temperature 

IfTiT in nuclear magnetic resonance[133, 134] (green diamonds). Lines are a guide to the eye. The phase transition temper¬ 
atures Ts . Tn and Tc[12] are indicated by thin colored lines and areas. , the temperature at which Cge extrapolates to 

zero, is indicated by blue triangles. 


experiment is very closely related to the bare nematic susceptibility Xv obtained using the shear-modulus 
data. In contrast, the values of Tq derived from the elastoresistivity data are considerably higher at all doping 
levels than those from the elastic and Raman data, which is unexpected within the Landau analysis of this 
quantity [49]. For this probe, Tq appears even to cross the Tg line at around 5% Co content, which is also not 
expected in the simple Landau theory. This curiously high value of Tq for the overdoped samples may partly 
be due to the exclusion of the low-temperature region from the fit of ruee- In this region, the divergence of 
rriQQ appears to be suppressed at lower temperature, which was attributed to disorder effects [128]. Such an 
effect appears to be either absent or much less pronounced in the Raman and Young’s-modulus data. The 
origin of the differences of the nematic susceptibility derived from the different experiments in both the Tq 
values for Ba(Fei_a;Coa;) 2 As 2 and detailed temperature dependence for optimally doped Bai_a;Ka;Fe 2 As 2 is 
unclear to us and deserves further attention. 


7. Magnetic correlations as the origin of nematicity? 

In the spin-nematic scenario of Refs. [3, 18, 19], magnetic correlations are at the origin of nematicity 
and, ultimately, the structural phase transition. Hence, it should be possible to derive the elastic properties 
of iron-based materials from their magnetic properties. Indeed, a scaling relation between 1/TiT, as a 
measure of the strength of spin fluctuations, and the shear modulus Cqq has been derived in this theory, 
which provides a useful test of the scenario [24]. Here, we briefly review the derivation of this scaling between 
spin-lattice relaxation rate l/Ti and Cqq and then we check the scaling using experimental data in three 
different iron-based systems. 
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7.1. Scaling relation between TiT and Cee 

An expression for the nematic susceptibility in terms of the dynamic spin susceptibility x is calculated 
in Ref. [19] as, 


(Xv) ^ 


1 


( 12 ) 


(see also Ref. [3]). Xv is renormalized by bilinear coupling to the elastic system (as in equation 1) to 


iXv) ^ 


1 

E,xHq) 


{go + A^/Cee.o), 


(13) 


where g = go + A^/Cee.o is the nematic coupling and the crucial parameter of the theory, q = (q,w) stands 
for the momentum and frequency dependence. The magnetic transition occurs when X^{g) diverges, but 

if (/ > 0 it is sufficient for Yl\x{g) to reach a finite threshold value (i.e., 1/g) to cause a divergence of x^ 
and induce the nematic/structural transition. This is the explanation why Tg can be higher than Tjq in this 
scenario even though both transitions are driven by magnetic fluctuations. The spin susceptibility y can be 
accessed by the spin-lattice relaxation rate divided by temperature, as measured in NMR experiments. 


1 


= Xa 


lim ^F 2 (q) 

CO —^ 


Imx(q,w) 


UJ 


(14) 


Here, loq is the NMR frequency, which is considered to be very small and F{q) is a momentum dependent 
form factor which peaks at the ordering wave vectors Qi and Q 2 when the magnetic field is applied par¬ 
allel to the ab plane [135]. As shown in detail in Ref. [24], IjTiT measured with the magnetic field in 
the ab plance is proportional to X)qX^('?) under certain approximations (assuming overdamped dynamics 
(x(q, w)) ^ = (x(q)) ^ — iojT, vicinity to a finite temperature critical point so that x'^ig) can be replaced 
by To X^(ci)) and the replacement of the form factor F(q) —> F{Q) because of the direction of the applied 
field). Using this proportionality, one can express Xv and, with the help of equation 3, Coe in terms of Ti, 


C'ee _ 1 

Coo.o ~ l + {aTiT-b)~^ 


(15) 


with the two parameters a and 5 [24]. The parameter b is particularly interesting, since it provides a measure 
of the nematic coupling strength b = -^^g. 


7.2. Test of the 1/TiT-Coo scaling relation in Co- and K-doped BaFe 2 As 2 

Figure 9a shows 1/TiT data of Ba(Fei_ 2 ;COa;) 2 As 2 from Refs. [133, 134], measured under an in¬ 
plane magnetic field and Fig. 10a shows the equivalent data for Bai_ 2 ,Ka;Fe 2 As 2 from Ref. [136]. We 
first discuss the Ba(Fei_a;Coa;) 2 As 2 system. Eq. 15 considers only the fluctuations around the AFM 
wave vector Q, which are referred to as “interband” contribution in Ref. [133]. The non-critical “in¬ 
traband” contribution has to be subtracted from the data as a background. Following Ref. [133], this is 
achieved by subtracting the data for the strongly overdoped Ba(Feo. 86 Coo,i 4 ) 2 As 2 , modeled as {1/TiT)^^^^^^ = 
0.11 K“^s“^ -1-0.63 K“^s“^ exp(—450 K/T) (black line in Fig. 9a). The obtained {^/TiT):^^^^^ can be scaled 
onto the elastic data according to equation 15 and the result is shown in Fig. 9b, with the scaling parameters 
a and b are given in the inset. The scaling works very well for all Co substitution levels, supporting a mag¬ 
netic origin of the shear-modulus softening[24]. Alternatively, a phenomenological linear scaling between Cqo 
and IfTiT was proposed in Ref. [137], which, however, does not work as well. Note that the parameter b is 
proportional to the difference between the Weiss temperatures Tg and Tq^TiT of the nematic susceptibility 
Xip and IfTiT, respectively, which are shown in Fig. 8 . Interestingly, Tg and Tq^TiT cross and b changes sign 
around the critical composition where the structural transition is suppressed. However, the assumption of a 
finite temperature critical point used to derive the scaling relation is not strictly valid in this doping region 
and, hence, this sign change of b might be an artefact. However, if b, or, equivalently, g really becomes 
negative for overdoped Ba(Fei_a;Co 2 ,) 2 As 2 , it would indicate that a non-nematic, C' 4 -symmetric magnetic 
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Figure 9: (a) Spin-lattice relaxation rate divided by temperature 1/TiT of Ba(Fei_a;Cocc) 2 As 2 with magnetic field H\\ab, from 
Refs. [133, 134]. The data for the 14% doped sample are taken as a background (black line), (b) Scaling analysis of the shear 
modulus Cqq/Cqq^o and spin-lattice relaxation rate. The lines are obtained from three-point bending, while the symbols show 
the data of (a), scaled according to eq. 15 after subtraction of the background. The scaling parameters are reported in the 
inset. 



Figure 10: Same as Fig. 9 but for Bai_a;KccFe 2 As 2 . NMR data with field H\\ab are taken from Ref. [136]. The vertical line in 
the inset marks the K content at which the structural transition is suppressed in a first-order-like fashion and the temperature 
dependence of the nematic susceptibility changes abruptly (see Fig. 6) 


19 


















state is preferred over the stripe-type one. Such a C' 4 -symmetric magnetic phase is, indeed, found to be 
induced by Na or K substitution in BaFe 2 As 2 very close to the critical point [ 86 , 92]. 

2 _ 2 

We note that a similar scaling attempt in Ref. [57], between the nematic charge susceptibility Xo ^ 
measured using electronic Raman scattering, and the elastic modulus failed. Here, the authors equated 

the Xo ^ with the renormalized nematic susceptibility x^ so that the scaling takes the form C'og/C'ee.o = 

( 2 2 \ — 1 2 2 
1 -I- axo ) (deduced using eq. 4) with a single parameter a. However, when assuming that Xo 

reflects the bare nematic susceptibility X(^[118, 131], the form of the scaling should rather be identical to eq. 
15, which does work well also for the Raman data of Ba(Fei_a;Co 2 ,) 2 As 2 . 

It is interesting to consider the scaling relation of 1/TiT and Cqq also in the Bai_a;K 2 ,Fe 2 As 2 system, for 
which both sets of data have been published[21, 136], but the scaling has not been attempted. The striking 
feature of this system is that both 1/TiT and Cge do not follow a Curie-Weiss temperature dependence 
over the whole doping region. In particular, the nematic susceptibility as obtained from the shear modulus 
changes its temperature dependence abruptly between 24% and 30% K content. The NMR data from Ref. 
[136] show a similar change of temperature dependence, however, occurring rather between 39% and 55%. In 
spite of the lack of a Curie-Weiss-like temperature dependence, it is remarkable that these data can be scaled 
quite well by Eq. 15 (Fig. 10). Tentatively, the same “intraband” background as for the Ba(Fei_a;Coa;) 2 As 2 
system has been subtracted from the raw 1/TiT data. The parameters a and b used in the scaling are 
shown in the inset of Fig. 10b. There are several points of interest here. First, there is an abrupt change of 
both a and 6 at ~ 25% K content (marked by a vertical line), reflecting that the temperature dependence 
of Cge but not of TiT changes abruptly, b is found to be close to zero just before the structural transition 
disappears, which is similar to the Ba(Fei_a;Coa;) 2 As 2 system. Beyond this concentration, b is negative 
indicating that the systems tends to a non-nematic magnetic order. Interestingly, the magnetic ground state 
of Bai_a;Ka;Fe 2 As 2 close to this K concentration seems, indeed, to be tetragonal and it would be fascinating 
to study NMR and shear modulus in detail for the respective substitution range. The parameter a also shows 
a strong dependence on K content, while it is roughly independent of Co content. Examining the scaling 
relation (Eq. 15) shows that the parameter a renormalizes the magnitude of TiT. The strong increase of a 
with doping reflects that the high-temperature values of 1/TiT increase significantly with K content, which 
might reflect either a change in the hyperflne coupling (as would be correctly captured by a [24]), or a 
strongly doping dependent background contribution. Note that l/TiT of the pure KFe 2 As 2 has the largest 
values of all samples [136], though it is supposedly far away from a magnetic instability. 

7.3. Shear-modulus softening and magnetic fluctuations in FeSe 

As described in the introduction, the iron-based superconductor FeSe is particularly interesting with 
respect to the relation of structure and magnetism. In particular, its large paramagnetic, orthorhombic (i.e., 
nematic) phase makes FeSe an interesting test case to study the origin of nematicity. As discernable from 
Figs. 4 and 5, Ijno] of FeSe is very similar to underdoped Ba(Fei_ 2 ,Coa;) 2 As 2 . In Fig. 11, we compare the 
two systems in detail. Curiously, hjno] and Xv of FeSe are nearly identical to Ba(Feo. 97 Coo,o 3 ) 2 As 2 , which 
has a similar rs[23]. The nearly identical temperature dependences indicate that the coupling A^/aCge.o has 
the same value in FeSe and underdoped Ba(Fei_a;Coa;) 2 As 2 , which is remarkable considering the differences 
between the two systems. 

In light of this similarity of the Young’s modulus, it is somewhat surprising to And that the evidence 
for spin fluctuations from the spin-lattice relaxation data of FeSe is much less pronounced than in the 
Ba(Fei_a;COa;) 2 As 2 system [23, 63] (Fig. 11b). In contrast to underdoped Ba(Fei_ 2 ;COa;) 2 As 2 , l/TiT of FeSe 
shows an increase upon cooling only below Tg and, actually, decreases upon cooling from room temperature 
down to Tg. This drastic difference already suggests that spin fluctuations may not be the origin of the 
diverging nematic susceptibility in FeSe. In order to exclude the possibility that a strongly temperature 
dependent “intraband” background contribution hampers the determination of the relevant (l/TiT);^^^^^. and, 
hence, the scaling analysis in FeSe, the NMR data have been analyzed in detail (see Fig. 12, [23]). 

Fig. 12(a),(b) show the spectral shift K°‘ and 1/TiT of a collection of ~ 10 single crystals of FeSe, 
measured in a held of 9 T with the held direction a both in-plane and along the c axis. K°‘ has a pronounced 
temperature dependence, which can presumably be explained in a Fermi-liquid picture by the small Fermi 
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Figure 11: (a) Young’s modulus Y[iio] FeSe, superimposed on the data for Ba(Fei_j;Coa;) 2 As 2 . FeSe fits very well in the 
series. The inset shows the nematic susceptibility (see eq. 3). (b) 1/TiT normalized by its (extrapolated) value at room 
temperature of FeSe and Ba(Fei_a;Coa;) 2 As 2 (data taken from Refs. [133, 134]). Vertical arrows indicate Ts. While 1/TiT 
of underdoped Ba(Fei_a;Cocc) 2 As 2 shows a strong increase on cooling towards Ts (and a steeper increase below), the data for 
FeSe only increase below T^. 



Figure 12: (a) NMR spectral shift and (b) spin-lattice relaxation rate divided by temperature \/T\T of single-crystalline 
FeSe with the field direction, a, in the ab plane (average of a and b axis in the orthorhombic state) and along the c axis, (c) 
shows a plot of ^ 1 IT\T vs. with the temperature as an implicit parameter indicated in units of K. The straight lines 

correspond to the Fermi-liquid type Korringa relation eq. 16, deviations from which show the emergence of a spin-fluctuation 
contribution. Such an additional contribution is evident only below Ts. The thin lines are a guide to the eye. 
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energy [138] found in FeSe [69, 70, 121, 122]. Since K and IjTiT are related by the Korringa relation for a 
Fermi liquid 

the temperature dependent K°‘ will lead to a temperature dependent 1/TiT. In order to isolate a possible 
contribution from spin-fluctuations to 1/TiT, we test eq. 16 by plotting y^l/TiT vs. with temperature 
as an implicit parameter in Fig. 12c. For both field directions, we find a linear relation between the two 
quantities, which shows that the Fermi-liquid Korringa relation is indeed satisfied down to Tg. Hence, 
this analysis shows the absence of any measurable spin fluctuation contribution to XjTxT for T > Tg, in 
strong contrast to underdoped BaFe 2 As 2 . This means that such spin-fluctuations cannot be the origin of 
the shear-modulus softening and, in consequence, the structural transition and nematicity in FeSe, and it is 
likely that the alternative orbital order drives the structural transition, as also suggested by ARPES studies 
[69, 71]. A similar conclusion was reached by another NMR work by Baek et al. [67]. Curiously, the nematic 
susceptibility Xv of underdoped Ba(Fei_a;Coa,) 2 As 2 and FeSe is very similar, which may raise doubt on the 
magnetic origin of nematicity in Ba(Fei_a;Coa;) 2 As 2 as well. Finally, the NMR data on FeSe (Figs. 11, 12) 
seem to suggest that the structural transition at Tg triggers the emergence of spin fluctuations. However, 
there is no such correlation under hydrostatic pressure, which enhances spin fluctuations [63] but suppresses 
Tg [76-78]. This result suggests an unusual relation between magnetic order and the structural transition 
in FeSe. 

As we pointed out previously [23], inelastic neutron scattering experiments are needed to determine the 
nature of the magnetic fluctuations, and very recently the first of such experiments have been reported 
on [72, 73]. Surprisingly, these studies provide evidence of magnetic stripe-like (7r,0) fluctuations in FeSe, 
similar to the 122 compounds. As in the other Fe-based materials, these magnetic fluctuations are found 
to occur already above Tg, in apparent contradiction to the NMR results. There is, however, evidence for a 
sizeable spin-gap in the neutron data already at 110 K, which can explain the absence of spin fluctuations 
in the NMR experiment, a probe that is sensitive only to the low-energy fluctuations. However, the neutron 
experiments, which have been performed both on polycrystalline [72] and single-crystalline material [73], 
are also not fully consistent with each other. In Ref. [72] the strength of the magnetic signal is temperature 
independent between 8 K and 104 K, whereas a strong temperature dependence is observed in Ref. [73]. 
These differences may also reflect differences in samples. The properties of FeSe samples depend strongly 
on the preparation technique. In particular, the samples prepared using a floating-zone technique used 
in Ref. [73] have undergone a phase transformation to the tetragonal, superconducting /3-FeSe phase on 
cooling to room temperature. Such samples are not single phase[139] and do not show a clear tetragonal- 
to-orthorhombic structural transition in the resistivity [139]. The samples prepared by solid-state synthesis 
used in Ref. [72] show an unusual difference of zero-field-cooled and field-cooled magnetic susceptibility 
around Tg\12], different from the vapor-grown crystals[64]. Hence, the issue seems not fully resolved yet and 
more careful neutron experiments on vapor-grown crystals are still desirable. Finally, we mention that these 
puzzling features of FeSe, i.e. the nematic phase and the absence of magnetic order, have attracted the 
attention of theorists, and there are now several different proposed theoretical scenarios, including strong 
magnetic frustration [82], spin-quadropolar order [81], the formation of a quantum paramagnet [83], and 
charge-current density wave [85]. 

8. Summary and outlook 

In this review we have compared the electronic nematic susceptibility of various iron-based superconduct¬ 
ing materials derived from measurements of the shear-modulus, the elastoresistivity, and the Raman response 
function. Particular emphasis has been put on our own studies of the Young modulus in Ba(Fei_a;COa;) 2 As 2 , 
Bai_a;Ka;Fe 2 As 2 and FeSe obtained via a three-point bending technique in a capacitance dilatometer. In a 
Landau formalism, in which an electronic nematic order parameter drives the structural transition via bilin¬ 
ear coupling to the orthorhombic lattice distortion, the nematic susceptibility was obtained from the elastic 
data. The relation of this “thermodynamic” nematic susceptibility to spin- and orbital degrees of freedom 
is found to be intricate, and may be different for the different iron-based systems. Notably, the nematic 
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susceptibility from the shear-modulus data seems to be closely related to the orbital nematic susceptibility 
from electronic Raman scattering in all the systems. On the other hand, the elastoresistivity behaves some¬ 
what differently, in particular for optimally doped Bai_a;Ka;Fe 2 As 2 . The excellent scaling of shear-modulus 
softening and spin-lattice relaxation rate 1/Ti in both Ba(Fei_a;Coa;) 2 As 2 and Bai_a;Ka;Fe 2 As 2 supports 
the notion that the structural transition is driven by magnetic fluctuations. FeSe seems to be an unusual 
case, in that there is no spin-fluctuation contribution to l/Ti above the structural transition, even though 
the nematic susceptibility of FeSe and underdoped Ba(Fei_a;Coa;) 2 As 2 have a very similar temperature de¬ 
pendence, suggesting the importance of orbital degrees of freedom. However, very recent inelastic neutron 
scattering experiments suggest that spin-fluctuations in FeSe are similar to the BaFe 2 As 2 -based systems. 
Further neutron studies on high-quality vapor grown FeSe crystals are thus highly desirable. 

The interrelationship of the various types of coupled order in iron-based systems—e.g., structural, orbital, 
magnetic and superconducting—has turned out to be a very rich held of study and electronic nematicity has 
grown into one of the most intensively studied concepts in the held. The relationship between nematic fluc¬ 
tuations and superconductivity remains an interesting and open problem. Recently, evidence for a nematic 
resonance in the superconducting state was presented in Raman scattering experiments [140] suggesting a 
close link between nematic fluctuations and superconductivity. Further, the enhancement of superconduc¬ 
tivity near a nematic quantum critical point has been investigated[141]. We believe that in order to obtain 
further insight into the microscopic origin of nematicity and its possible relation to superconducitivity, the 
study of systems in which structural and magnetic transitions do not closely follow each other is particularly 
promising. In this context, FeSe and the magnetic 6*4 phase in hole-doped BaFe 2 As 2 have recently attracted 
a great deal of attention and may very well hold further surprises in the near future. 
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